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Abstract, Individual particle trajectories are solved for 


in three models of possible field configurations of the geomagnetic 

tail- The analytical results of part 1 [ Speiser , 1965 b] are 

applied to two models. Both models contain magnetic field lines 

oppositely directed on either side of a neutral sheet, with an 

electric field perpendicular to the magnetic field and parallel 

to the neutral sheet. The models differ in the rate of variation 

of a magnetic field component perpendicular to the neutral sheet, 

and hence in the rate of field -line-crossing of the neutral sheet. 

For the two models, particles are accelerated and turned toward 

the earth within the neutral sheet and are ejected from the neutral 

sheet with small pitch angles to a magnetic line of force, with 

energies of tens-of-kilovolts . For the first model, a dipole-plus- 

tail model, electrons are ejected at about 150 and protons 

about 50 R back in the tail. For the second model, an extended - 

tail-model, electrons are ejected at about 500 R , and protons at 

about 400 R . Proton auroras would be expected about l/2° lower 

latitude than electron auroras, and isotropic fluxes should be 

measurable out to distances of the order of 2.5 R from the earth. 

e 

Extremely thin sheets of incoming particles are produced, about 
1 km for electrons between 1 and 10 kev. The third model consists 
of a three-dimensional dipole field added to a tail/neutral -sheet 
field, and trajectories are calculated numerically, finding spatial 
regions of high intensity using Liouville 1 s theorem. These 
spatially intense regions are near the auroral zones when mapped 



onto the earth; they move to lower latitudes on the earth with 
an increase in the strength of the tail field, and their thicknes 


is roughly proportional to the thermal velocity of the particles 
incident on the tail. The models may be applicable to other 
situations where neutral points or sheets may exist, such as the 
day-side magnetospheric current sheet, the interplanetary field. 


solar flares, etc. 
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Introduction . This paper is concerned with the possibility 
that the geomagnetic tail and more specifically, the neutral 
sheet in the tail, accelerates particles from the solar wind and 
is the immediate source of auroral particles. 

Parker [1957] has investigated particle motion about a 
neutral sheet, and has found that regardless of the initial particle 
configuration, stability soon results with the current given by 
curl B, just as in classical hydromagnetics. Weiss and Wild [1964], 
and Chapman and Kendall [1964] have also looked at particle trajectories 
about a neutral line without assuming the existence of any electric 
fields . 

An essential part of Dungey’ s theory, however, concerns an 
electric field. He showed [ Dungey, 1958] that if one assumes 
infinite conductivity at a neutral point, an infinite discharge 
develops in a finite time. Any finite conductivity would therefore 
limit the discharge and an electric field across the neutral 
region would remain. Such an electric field would be consistent 
with the "frozen-in" field of the advancing plasma. 

Akasofu and Chapman [1961] developed a neutral line discharge 
theory of auroras, which is based upon one or more ring currents 
and with which they also hoped to explain the main phase of magnetic 
storms. Although this work was specifically pointed toward the 
ring current formation of neutral lines, later they suggest 
[ Akasofu and Chapman , 1964] that the ring current is just one of 
the qualified sources of neutral lines and that in the general 
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sense, "the field line projection of an arc in the ionosphere on 

the equatorial plane must be a 1 singular 1 line in the sense that 

on the line the guiding center approximation breaks down." 

Piddington [i 960 ] suggested that if some field lines from 

the earth were trapped into the solar wind at the time of a 

geomagnetic storm, a geomagnetic tail would be formed. His 

suggestion was made primarily to explain the main phase of geomagnetic 

storms, but he also mentioned that it might have "incidental effects" 

including auroral and Van Allen radiation, the gegenschein, and 

the diurnal variation of cosmic rays. 

The neutral sheet in the geomagnetic tail has been probed with 

the Imp-1 satellite [ Ness , 1965 ]. The sheet is about 0.1 

(Earth-radii) thick and extends from about 10 R to at least 30 R 

e e 

(the apogee of Imp-l), and is found on almost every pass through 

the tail. Cahill 1 s [1964] results using data from Explorer XIV 

complement Ness 1 results and show the tail structure in closer to 

the earth. These results generally fit the pattern of dipole field 

lines being stretched out in the tail and reversing across a 

neutral sheet beginning somewhere around 10 R . 

e 

McDiarmid and Burrows [ 1965 ] using data from the Alouette 1 
satellite at 1000 km, show occasional spikes of high intensity 
electrons (energies above 40 kev, fluxes approaching 10 9 cm^ 2 sec'‘ 1 ster" 1 ) 
occurring in a narrow latitudinal range at latitudes above the 
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trapping region. These field lines are therefore presumably 
connected to the tail and the tail is the source of the particles. 
Many other particle measurements in the tail show regions, usually 
near the neutral sheet, where the electron temperature is high, 
particle densities are much higher than interplanetary values, 
and intense fluxes or "islands" of energetic electrons are seen. 

[See Anderson, et al. , 196 5 ; Anderson, 1965; Frank, 1965; Gringauz, 
et al., I96O; Howard, et al . , 1965]. 

If the magnetic field, B, at the center of the neutral sheet 
goes to zero, the gyro-radius increases without bound, so no 
matter how slowly B may vary, R/L (L, a characteristic system length) 
will be large, and adiabatic theory cannot be used. Particle 
trajectories must therefore be either solved for analytically from 
the equations of motion or computed numerically. If B is small 
but non-zero at the center of the sheet, R/L must be determined to 
see if adiabatic theory can be used. 

In the following section, particle trajectories are calculated 
applying the results of the non-adiabatic analytic theory [ Speiser, 
Part 1 , 1965], (hereafter referred to as Si), to two models of 
the fields expected in the geomagnetic tail. A section follows 
where numerical trajectory results are given for a more complicated 


magnetic field geometry. 
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Analytical Applications . Although the geomagnetic tail has 
been probed, and magnetic field measurements have been made by 
Cahill [1964] out to about 12 R^ and by Ness [1965] out to about 
30 R^, measurements have not been made further back, and the exact 
description of the magnetic field inside the neutral sheet is not 
known. (An exception is the null result of Mariner IV [ Coleman et al. , 
1965] where the geomagnetic tail was not found at 3300 R^.) A simple 
model of the fields in the geomagnetic tail will be assumed as 
indicated in Figure 1. 

The electric field is taken to be uniform and across the tail 
(in the -e^ direction of Figure l). Such an electric field is based 
upon the flow patterns of the open model, [ Dungey , 1961; 1958] and the 
work of Axf ord and Hines [1961], Axford , Petschek and Siscoe [1965], 
and Petschek [196k], and on ground observations of high latitude 
magnetic field fluctuations and their inferred current and equipotential 
systems. If an electric field does exist across the tail. It would 
most certainly have time and spatial dependencies. Such variations 
are not included in this model, but would modulate the results 
presented here. Measurements of electric fields In the magnetospheric 
tail have yet to be made, hence questions as to their existence are 
argumentative . 

Using the fields of Figure 1, the results of SI can be used. 

Those results are summarized as follows. 
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1. Particles of either sign incident on the neutral sheet 
oscillate about the sheet due to the reversal of the magnetic field. 

2. As a particle oscillates about the sheet, it gains energy 
from the electric field, and is turned into the -e direction (toward 

y 

the earth) by B^. 

3* A particle will oscillate until it has been turned so much by 
B that z, its velocity in the e direction, changes sign. At that 
time the particle is ejected from the neutral sheet. (The ejection 
time is t = n/( q/m)bT]. ) 

4. The velocity at ejection is almost entirely in the -e direction 
if T] « 1. (The ejection velocity y(T) = -2a/l|b.) Thus the ejection 
pitch angle (a) will be small if T] is small, (a Tj/2 jx 0 /(a/b) - 2j . ) 

5 . In the moving system where the electric field is transformed 
away, the motion is seen as an oscillation in x about the neutral 
sheet combined with a circular drifting of the trajectory in the 
neutral sheet (the y-z plane). The neutral sheet effectively uncouples 
the circular drift from the oscillation about the sheet. (Although 
the oscillation is. coupled to the circular drift.) 

The electric field strength in the tail, "a", is assumed to be 
about 0.3 volts/km. This gives a potential of about JO kilovolts 
across a tail of 40 R diameter. This seems to be the right order 
of magnitude for the potential across the polar cap, which is mapped 
into the tail, assuming magnetic lines of force are e qui potent ial s . 

Since particles gain energy by drifting across the tail while 
oscillating about the neutral sheet (item 2 above), an electric 
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field of different magnitude will affect the ejection velocity 
(item 4 above) , ana will limit the maximum attainable particle 

energy . 

The magnetic field strength, "b", is taken to be 20Y 
(ly = IQ” 5 gauss = 10“ 9 weber/m 2 ). Ness [1965 ] finds the 

field to be from 10 to and about 40y at the time of a magnetic 

storm. 

The other parameter, T], that we need to know is the ratio of 

the magnetic field component perpendicular to the neutral sheet, B , 

x 

to "b", the solar-antisolar field strength outside of the neutral 

sheet. As a first case, we will assume that B is furnished by the 

earth’s dipole field as sketched schematically in Figure 1. This 

assumption is certainly artificial and merely provides a model for 

the rate of crossing of field lines through the neutral sheet. Using 

such a model implies that the field strength at the center of the 

neutral sheet at 16 R^ is about yy 3 and this seems to be larger than 

that measured by Ness [p. 2993* Figure J, 1965 ]* Another model for 

B is discussed at the end of this section, 
x 

Referring to Figure 1 for the values of the constants, and using 
the results of SI, Table 1 can be constructed. Figure 2 shows a 
sketch of particle trajectories for this model. 
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The analytical study (Si) was based on the assumption that Tj is 
a constant. From Table 1, for protons at 50 R e > it is seen that 
the particle drifts 25 R^ toward the earth before it is ejected from 
the sheet. Using the dipole model, 7] would change by a factor of 
eight, so the above-mentioned approximation does not seem very good. 

The larger T| would serve to turn the proton in a tighter circle, 
and the proton would be ejected sooner and would gain less energy 
than indicated in the table. However, if a trajectory is broken into 
a number of small segments, over each of which 7] is approximately 
a constant, then SI can be used over each segment. For the above 
example, when the proton has drifted in to about 40 R^, it has 
already gained about 20 kev, so the inconstancy of 7] does not affect 
the results as much as would be at first supposed. This is because 
the biggest part of the energy gain is during the first part of the 
trajectory before the particle has been turned much in toward the 
earth. 

In SI the qualitative behavior of the oscillation about the 
neutral sheet has been determined, but the details have not been 
solved for analytically. Knowledge of the output pitch angle requires 
the detailed knowledge of all of the velocity components at output, 
x was estimated as of the order of x Q , and z was estimated as of the 
order of zero, since ejection occurs when z changes sign. The pitch 
angles in Table 1 are therefore shown for several values of x Q , making 


the above assumptions. 
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Many proton trajectories have been computed, solving the 
equations of motion using a Runge-Kutta computer subroutine. Using 
"j - l/by (from the earth's dipole field at 50 R g ) and the values of 
"a" and "b" as used in Table 1, it is found that ejection pitch 
angles lie between 0 and 6 degrees for incoming velocities of zero to 
200 km/sec. 

From Table 1, at 50 R g if x 0 !=a 12u « 200 km/sec, a pitch angle 
of about 4° is expected, and is therefore in general agreement with 
the computed value. However, the computations do not show the pitch 
angle to be related to x G in as simple a manner as found from the 
qualitative arguments for determining the pitch angle of SI (paragraph 
4. above). 

The analysis (Si) of the energy gained, the turning of the 
trajectory toward the earth, the trapping in the neutral sheet and 
subsequent ejection when z changes sign are all confirmed by the 
trajectory computations. The qualitative behavior, oscillation frequency, 
amplitude variation, etc., about the neutral sheet also agrees with 
the theory. At ejection, however, the perpendicular velocity components 
x and z, and hence the pitch angle, depends on where the particle 
is and what its velocity is when z changes sign. If the particle is 
close to the peak of its last oscillation when z changes sign, it 
will have a small x and thus small pitch angle; if it is close to 
x = 0, its x will be large and will increase until ejection, and will 
thus imply a large pitch angle. 
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Figure 3 is a machine computation of a proton trajectory 

using the fields of Figure 1 and the constants as in Table 1. The 

initial velocities at the neutral sheet (x = 600 km) are indicated. 

The comparisons of the theoretical (Si) predictions with the computed 

values are indicated on Figure 3* 

Dungey [1965] predicts a tail length of the order of 1000 R^, so 

a neutral point may exist at about 500 R^. A possible model for the 

perpendicular component, B , would therefore be one which goes 

x 

linearly from about 1Y at about 50 to zero at 500 R & . This 
model is sketched in Figure 4. The limiting field line from the 
neutral point as sketched is attached to the auroral zone in agreement 
with Dungey 1 s [1961] open model. Table 1 can be used for the 
application of the analytical results to this model, with the only 
modification being the first entry, that is the distance back in the 
tail at which the particle enters the neutral sheet. For the first 
column B is l/4Y at y = 388 R , so protons of about 30 kev will 

X 6 

come from this region, and electrons of about 12 kev energy will 
come from the region where B^ 10" 2 V or from about 495 R e * 

From conservation of flux we can find the latitudinal separation 
at the earth of the two field lines which come from 388 R^ and 495 R 0 - 


That is: 
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r B-dA = r B-d A , or (l) 

earth tail 


01 ^95 R 

B 0 R 0 2 f sine d 0 = f B x (y) T d y (2) 

23° 386 R g 


where B 0 « 60, 000Y (assumed constant if d l is small), B (y) is 
taken from Figure 4, and T, the tail width, is assumed about 40 R . 
For the above example, 9 X comes out to be about 23*4°; thus for 
this model of the fields, protons are ejected along field lines which 
intersect the earth 1 s surface about 0.4° lower latitude than do the 
field lines along which energetic electrons are ejected. 

From the analytic study, SI, it was found that the velocity of 
the particles ejected from the neutral sheet varies inversely as 
B^. 2 . For electrons, the relationship is therefore 


W = (1.2/B 2 ) 

X 
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where W is the energy in ev, and is in gammas. For the model 
used in Figure 4, we have 


B 

x 


5QQ- y 

450 


Y 


( 3 ) 


with y measured in earth-radii (R ). Therefore the ejected electron 

e 

energy as a function of distance back in the tail is 


2.4 x 10 5 eV 
(500-y) 2 


(4) 


and the equivalent proton energy is found by multiplying the right- 
hand side of equation (4) by 1 , 836 . Both expressions for the energy 
are valid until the maximum potential across the tail has been gained 
(see Table l) . The co-latitude is found as a function of distance 
in the tail by conserving flux as before, and from equation (4), 
the electron energy can be found as a function of co-latitude, and 
this is : 


0.14 


9 - 9o 


w 


(5) 
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where the angles are in radians for W in eV, and 0 O = 25°. From 
equation (5) the latitudinal separation on the earth for ejected 
electrons from 1 kev to 10 kev, for example , is found to be about 
1.5 x 10” 4 radians, which corresponds to a beam width of about 
0.8 km thickness at the earth. 

M. P. Nakada [private communication, 1965] has suggested that 
field line loading may be an important problem for any auroral theory 
in which the particle source is located far away. That is, fluxes of 
particles with energies and intensities large enough for auroras 
may have more energy density than magnetic field energy density if 
they come from weak regions of magnetic field, and if these fluxes 
are isotropic in the weak -field-region. 

O' Brien [196M observed that the average energy flux of electrons 
with energies greater than 1 kev precipitated In the auroral zone is 
about k ergs cm“ 2 sec" 1 . He also notes that fluxes may occasionally 
be as high as 2000 ergs cm~ 2 sec” 1 , [ Mcllwain , i960; 0 1 Brien and 
Laughlin , 1962]. Assuming these fluxes are of 10 kev electrons, the 
above numbers imply fluxes of 3 x 10 8 to I.5 x 10 11 electrons cm" 2 sec" 1 , 
and thus densities of .05 to 25 electrons cm" 3 . Omholt [1965] 
estimates densities as high as 1000 electrons cm" 3 at the time of an 
intense aurora. 

To see what sort of densities this implies for a source, assume 
for the region between the source and the earth Liouvilie 1 s theorem 
holds; energy is conserved (after ejection for the model Si); p, the 
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first adiabatic invariant is conserved. Liouville' s theorem allows 
us to write 


J (particles cm" 2 sec -1 ster"* 1 energy"* 1 ) 
E 


constant (6) 


which implies J = constant from the second assumption above. The 
number density N, (particles cm" 3 ) is found by integrating Jv^ 1 
(v is the particle velocity) over solid angle and energy. In this 
example it will be assumed that J is independent of the solid angle 
Q, (the flux is isotropic) , and exists for only one energy (the 
flux is monochromatic). Then, at the earth, 


N 

e 


J 6(E -E 0 ) dQ dE 
e e e e 

( 2^)1 

m 


2E^i 


(£££) 

m 


( 7 ) 


where the subscript e refers to the earth. In the tail (subscript t) 
we have 


N t = J 


J t 6(E t - Eo ) dn t dE t 




■» 


dn t 

W 


dQ 


2E, 

in 


i 


( 8 ) 
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where use has been made of Liouville' s theorem. To evaluate N^, we 
use the conservation of u and E to obtain: 


sin 2 or e 

”b 

e 


sin 2 ot. 


B, 


( 9 ) 


along with the definition of solid angle 


dQ^ = 2tt sin d 

( 10 ) 


dfi = 2tt sin o- d a 
e e e 


where o' is the particle pitch angle. Therefore, 


dQ t 

dQ 

e 



B 

e 


cos a 


e 


J 1 - B. B _i sin 2 -a 
v t e e 


( 11 ) 
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and, evaluating the integral. 


2ttJ 

N, = ~t (1 - ^V^B7B~ =r ) ( 12 ) 

t (2E 0 /m)2 t e 


so 

N,_ 

(13) 
e 


IT = ( 1 - ■J 1 - V B e > 


or 


N t 

tT = h B, /B for B « B . 
N t 7 e t e 


(no 


If we now equate particle energy density in the tail to magnetic 
energy density, we have 


E 0 N t = 2.65 B t 2 


( 15 ) 


for Eq in eV, and B in gammas, and using equation (l4) we have 
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B 


t 


Ep 

5-3 B e 


(16) 


or, for E 0 - 10 kev, and B g = 60, 000y, 


B, 


Ke 

32 


( 17 ) 


Thus to have the same energy density in the field as in the particles 
for O'Brien's 2000 ergs cm" 2 sec -1 (N e = 25 electrons cm -3 ) requires 
a tail field of only about ly. If Omholt' s estimate of « 1000 
electrons cm -3 is at times correct, a tail field of about 30 Y would 
have about the same energy density as particle energy density. If 
we make the further assumption that the particle energy density be 
at most 10$ of the tail magnetic energy density, then tail fields of 
10 to 300y are required. If, however, auroral particles are shot 
nearly along magnetic lines of force from the source, as they are from 
the mechanism reported in this paper, then the pressure perpendicular 
to magnetic lines of force, tending to push field lines out can be 
much smaller than the isotropic case. Therefore a smaller B^ can 


contain the particles . 
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Discussion of Analytical Applications . The significant results 
of the previous applications of the analytic theory (Si) to a tail 
model of the earth are that particles are accelerated, turned in 
toward the earth and ejected from the neutral sheet with small pitch 
angles to the magnetic lines of force. This ejection at small pitch 
angles may be important for the prevention of "field line loading" as 
discussed in the previous section. Certainly more experimental 
observations on particle fluxes and energies and magnetic field 
configurations are needed before some of these questions can be 
answered. 

These applications also predict proton auroras to be at lower 
latitudes than electron auroras, (for the second model the latitudinal 
separation is about 0.4°). Omholt [1963] says "There is often 
(perhaps always) a distinct dark region (up to 1° in latitude) 
between the ’proton aurora' and the main forms." 

Electrons in the energy range from 1 to 10 kev would be found in 
an extremely thin beam at the earth, i.e. about one kilometer 
thickness. Equation (5) also predicts a hardening of the beam with 
latitude, but 'the electric drift between the ejection point and the 
earth has been neglected, (this would tend to move the higher 
energies to slightly lower latitudes), and the self-consistency of 
such a thin beam has not been investigated. This result may therefore 


be incorrect. 
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The analytic results also predict a monochromatic beam at a 

given point in space, but those results are approximate and based 

on the assumption that the perpendicular component of magnetic 

field in the neutral sheet, B , is constant. For a linear variation 

of B as used here, a different spectrum should be found, but this 
x 

requires computing, and has not yet been done. 

0 1 Brien [1964] suggests that a major experimental study should be 

made to determine the limits of isotropy in the incident beam and 

thus the cause of auroral precipitation. He finds that for electrons 

with energies greater than 40 kev, fluxes are isotropic to within 

10$ at 1000 km. The results of the present study imply fluxes from 

the tail within a cone of the order of 0.1 radian, and thus one would 

expect to measure isotropic fluxes out to about 2.5 earth-radii. 

The difference between the second model used, (as in Figure 4), 

and the dipole-plus -tail model lies in the distance from the earth 

at which particles are ejected from the neutral sheet. Indeed, any 

model giving a different variation of B (hence a different rate of 

X 

crossing of field lines through the neutral sheet) -will merely move 

the particle ejection points toward or away from the earth. If, 

however, B is negative, then particles will be turned away from 
x 

the earth and will be ejected into the tail (in the anti-solar 
direction) . Close to the earth, however, B^ is positive, so it is 
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likely that it becomes very weak far from the earth, and it would 
reverse if there is a neutral point somewhere in the tail. Such a 
reversal is likely if the interplanetary field has a southward 
component and the field has to eventually fit onto the interplanetary 
field. Thus, if there is a neutral point in the tail as indicated 
in Figure 4 , and if the electric field is as indicated, particles 
will be shot in toward the earth between the earth and the neutral 
point, and will be shot away from the earth on the other side of 
the neutral point. 

If the electric field in the tail is radial as Taylor and Hones 
[1965] suggest, then this mechanism can work over the dusk half of 
the tail (meaning negative values of z), but will not work over the 
dawn half (positive z) . Their electric field is in the 4 -e^ direction 
for the dawn half of the tail, and particles will not be trapped in 
the sheet in this region. However, the bulk flow velocity or 
circulation is such as to bring particles into the dusk half for their 
model . 

For particles to be trapped in the neutral sheet, it is only 
necessary that the magnetic field reverse across the sheet so that the 
magnetic force term in the lorentz-force equation reverses across the 
sheet. The specific linear variation used in SI lends itself to 
analytic solution and is probably valid over some portion of the 


;heet . 
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The self-consistent problem has not been solved. However, one 
interesting result of SI and these applications is that particles 
of different masses injected at the same point on the neutral sheet 
are trapped, accelerated, turned, and all are ejected with the same 
velocity. So from this simple viewpoint, space charges should not 
build up to invalidate the mechanism. This, however, assumes a 
neutral sheet of infinite extent, and is probably correct over the 
regions where the protons are turned and ejected toward the earth with 
energies in the tens of kilovolts. But in the region where the electrons 
are ejected with the same energies, protons drift completely across 
the tail, and thus leave the sheet (when the maximum potential has 
been gained) with smaller velocities than the electrons. Space 
charges may thus build up in this regime, limiting the electric field. 
Whether or not this is a real problem must await a self-consistent 


analysis . 
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Numerical Solutions . Speiser [1965a], presented results of 
numerical solutions of particle trajectories about a model current 
sheet. The model used was chosen for the fields "guessed" in the 
geomagnetic tail, but no connection was made to field lines at the 
earth. Those results showed that accelerated protons, emergent along 
magnetic lines of force, have their greatest intensity in a thin output 
sheet. A similar study was not made for electrons because of the 
enormous amount of computer time involved, but on the basis of a few 
trajectories it was suggested that the spatial structure for electrons 
would be much thinner than that for protons. These results are 
substantiated by the analytical results in the previous section of 
this paper. 

In order to make mappings of the output sheets onto the earth, 
a dipoie-plus-tail model is used. This model is the same as that 
used for the analytical study in the previous section, but ail of the 
components of the dipole magnetic field are kept, and not just B as 

X 

in Figure 1. The method of procedure is the same as used before 
[ Speiser, 1965a] and will only be summarized here. 

1. A three dimensional dipole field is added to 
a tail field which is as indicated in Figure 1. This field is 
certainly incorrect in the dayside magnetosphere and will only be 
used for calculations in the tail. 

2. A proton trajectory j_ s started at some point on 
an output plane (x 0 -z 0 plane), with a velocity directly along the 
magnetic line of force through that point, and the trajectory is then 
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numerically solved backwards in time. After the particle has passed 
through the neutral sheet (backwards in time) its velocity on the 
imput side is noted. Assuming a Maxwellian distribution, a value for 
the distribution function can then be found on the input side. 

5- Using Liouville' s theorem, the distribution 
function is constant along a trajectory; an intensity map over the 
output plane (xo-zo) may then be made. 

4. Regions of greatest intensity may then be 
mapped onto the earth by solving for the trajectories forward in 
time from the output plane. 

Figure 5 shows backwards plots for three protons, with only the 
x-component versus time. The strength of the tail field, b, is 20y, 
the y and z starting positions are Yo = 26.5 R g , z Q - 0.1 , 

and v 0 = 1664 km/sec for all trajectories. It is seen that if the 
particle starts with too large an x Q , it follows a field line back for 
a time, and then the electric drift causes the particle to drift away 
from the neutral sheet (backwards in time). Since the particle's 
energy changes only in the neutral sheet, its speed is the same on the 
input side (t s* -170 sec). A second particle starts at Xo « 3000 km, 
reaches the neutral sheet and oscillates about 16 times before input. 

Its energy has thus changed very much and the particle can have a small 
velocity at input. The remaining trajectory also reaches the neutral 
sheet, but does not stay in long before input (t « -50 sec) . Its energy 
has therefore not changed much, and its input velocity is not as small 
as the preceeding particle. Using the ideas from SI, we can say that 
this latter particle did not stay so long in the neutral sheet because 
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was larger for it and the particle was ejected sooner than the 
preceeding particle. 

If one therefore expects the majority of particles incident on 
the neutral sheet to have small velocities, then the middle trajectory 
of Figure 5 should indicate the region of the largest output intensity, 
since particles with nearby trajectories come from highly populated 
regions of velocity space. This is a crude explanation of the use of 
Liouvilie' s theorem, and of why such a model produces the largest 
intensities of output particles in thin sheets. 

Writing the distribution function as 



where F = v - u, v is the particle velocity, u is the bulk flow velocity, 
and v_j. is the thermal velocity, we can see that f is largest when 
F < v . 

Figure 6 shows a contour map of F on the x 0 -z 0 plane. F does 
reach lows of about n 100 km/sec" in regions near the center of the 
M 200 km/sec M contours, and the breaks in the n 200 km/sec" contours 
are probably not real, i.e., they would probably disappear if a plot 


with a finer net were used. 
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Figure 7 shows the positions of the output sheets with large 
intensity (F < v , and « 100 km/sec) on the x 0 -z 0 plane for 
various tail fields. 

Figure 8 maps the intersections of Figure 7 forward onto the 
earth, using the dipole -plus -tail model. 
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Discussion of Numerical Results . These numerical calculations 
have only been done for protons using a three dimensional dipole - 
plus -tail model. The results show that the largest intensities of 
particles accelerated in the neutral sheet and coming out along 
magnetic lines of force occur in thin sheets. (If v^ is taken to be 
larger than 100 km/sec, from Figure 6, the thickness of the output 
sheets is increased.) These sheets when mapped onto the earth 
intersect the earth near the auroral zones, and the sheets move 
to lower latitudes with an increasing tail field. These results 
are for only one output velocity, 1664 km/sec corresponding to a 
proton energy of about 12 kev. A longitudinal assymmetry would 
be expected for higher energy protons, in a plot like Figure 8, 
since the tail has only finite width, and the particles must drift 
across the tail in the neutral sheet to gain energy. 

Figure 8 is not expected to be quantitatively accurate since 
the dipole-plus-tail model is not very good near the earth, and 
no good on the day-side. Williams and Mead [1965, page 3025, Fig. 6] 
develop a magnetic field model in the magnetosphere which utilizes 
the dipole field, the field due to magnetospheric surface currents, 
and a current sheet in the tail. Using their Figure 6, the field 
line passing through x 0 = 3000 km « l/2 R g , at y Q “ 25 R g comes 
from a latitude between 65 0 and 70°. Therefore, using their model 
for mapping the intersections of Figure 7 onto the earth would lower 
the mappings of Figure 8 by about 5°> which would be in better 
agreement with the auroral zones. 
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Summary and Conclusions . A simple model is used for the 
possible electric and magnetic field configuration in the tail of 
the earth's magnetosphere. Satellite measurements out to about 
30 have shown the general result of high latitude field lines 
being pulled back in the tail, with the formation of a neutral-or 
current-sheet across which the field reverses direction. Not much 
is known yet about the field at the center of this neutral sheet, 
or about the existence of any electric fields. For the models 
discussed here, an electric field across the tail generally from 
dawn to dusk is assumed to exist, and a weak magnetic field component 
perpendicular to the neutral sheet and tending toward zero with 
distance is assumed inside the neutral sheet. These features need 
not be permanent, but can cause particle trapping and acceleration 
in the neutral sheet whenever they may exist. The self-consistency 
of these solutions and the effects of various instabilities should 
be investigated but this has not yet been done. 

Two models, a dipole -plus -tail model (Figure l) and an extended 
tail model (Figure k) are used to apply the analytic theory (Si). 

In both cases, electrons and protons incident on the neutral sheet 
are trapped, accelerated, turned toward the earth, and ejected frcm 
the neutral sheet with small pitch angles to a magnetic line of force. 
The main difference between the two models is the ejection point 
back in the tail. For the two models electrons of the order of 
10 kev are ejected at about 150 R g and 500 R g respectively, and for 
protons of about 30 kev, the ejection distances are about 30 R e and 
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400 R , respectively. The model produces extremely thin beams of 
incoming particles at the earth, about 1 km for electrons between 
1 an d 10 kev. Proton auroras are produced at latitudes about l/2° 
lower than electrons, and isotropic fluxes should be expected out 
to about 2.5 R • 

Many trajectories of about 12 kev protons have been computed 
for a three-dimensional dipole -plus -tail model. Using Liouville' s 
theorem, regions of large intensity are found which are near the 
auroral zones when mapped onto the earth. An increase of the tail 
field moves the mapping onto the earth to lower latitudes. 

These features agree generally with some auroral observations, 
namely: the latitude of the auroral zones; the energies of auroral 

protons and electrons; the fluxes of auroral protons and electrons; 
the appearance of proton auroras at lower latitudes than electron 
auroras; a dawn/dusk assymmetry for electron/proton auroras; the 
movement of auroral forms to lower latitudes with an increasing 
tail field (hence solar wind pressure); the thinness of auroral 
forms; and the gross conjugacy of auroral events. 

The models may be applicable to other situations where neutral- 
points or sheets are thought to exist, such as the day-side 
magnetospheric current sheet, neutral points in the interplanetary 
field, the influence o'f Jupiter' s satellite Io on Jupiter' s radio 
emission [ Warwick and Dulk, 1965 ], and solar flares. Since the 
models predict particles accelerated and ejected with the same 
velocity independent of mass, a comparison could be made, for example. 
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of helium and proton energies in a flare (or in the geomagnetic 
tail) to see if they are in the ratio of the masses. If some of 
the recently observed "red spots" on the moon [ Kozyrev, 1959 ; 19 & 3 ; 
Greenacre , 1963; Cameron, 1965] are due to fluorescence, the particles 
could come from this mechanism while the moon is in or near the 
earth' s tail. 
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FIGURE CAPTIONS 


Figure 1. Dipole -plus -tail model (in the meridian plane containing 
the earth-sun line). The thickness of the neutral sheet is 2d, 
the magnetic field strength outside the sheet is "b", 

or "T]b" is determined from the dipole field at a particular 
location, and the electric field strength is "a". Arrows 
are magnetic field components. 

Figure 2. Sketch of particle trajectories using the fields of 
Figure 1. Both protons and electrons oscillate about the 
sheet accelerating in opposite directions, and are turned in 
the same direction by the small magnetic field component 
perpendicular to the sheet. The direction turned is toward 
the earth if this perpendicular component is northward 
as indicated, otherwise into the tail if the component is 
southward. When the particles are turned 90 ° , they are 
ejected from the neutral sheet. Electrons come out much 
sooner than protons, with the same velocity as protons, 
hence gain less energy. The dimensions shown are illustrative 
and not to scale. 

Figure 3. Isometric trajectory plot of a proton trajectory in 

the dipole -plus -tail model, see Figure 1, B = lAv, b = 20y, 
d = 600 km, a = l/3 volt/km. For this model the earth would 
be at -50 R . Initial conditions; x b 600 km, y = z = 0, 


T is the 


x = -60 km/sec, y = 15 km/sec, z = 10 km/sec. 
ejection time, a is the ejection pitch angle. 

Figure k. Extended tail model. Field lines are sketched for the 

above model, where B is assumed to vary as indicated, going 

x 

from ly at 50 R to zero at 5 00 R . 

e e 

Figure 5 • Backwards trajectory plots, starting at z 0 = 0.1 R , 
yo = 26.5 R , v Q = 1664 km/sec. 

Figure 6. F contours in units of 200 km/sec, y 0 = 26.5 R • 

Figure 7* Intersections of sheets with largest intensity with the 
Xo-Zo plane (y Q = 26.5 R ) . 

Figure 8. Mapping of intersections of Figure 7 onto the earth 

using the simple dipole -plus -tail model. These mappings move 
about 5 ° lower in latitude when Williams and Mead 1 s [1965] 
field model is used. An asymmetry would be expected for 
higher energy particles, or for the same energy if the electric 
field were weaker, when a tail of finite width is considered. 
Protons would be shifted toward dusk and afternoon, and 
electrons would be shifted toward dawn and morning. 
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